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The effective pair potential between mesoscopic cliarged particles in a neutralising background 
medium takes a Yukawa form exp(— Ar)/r with screening length A~^. We consider a dilute suspen- 
sion of such Yukawa particles dispersed in a solvent with correlation length ^ < and show that 
the Yukawa interaction is 'screened' if the pair potentials between solvent particles exhibit Yukawa 
decay with the same screening length A"'^. However, if the solvent pair potentials are shorter ranged 
than the solute Yukawa potentials, then the effective potential between pairs of solute particles is 
'unscreened', i.e. the effective potential between the solute particles is equal to the bare potential 
at large particle separations. 



It is well known that the effective potential between 
charged particles immersed in a neutralising mobile 
medium is shorter ranged than the 'bare' Coulomb inter- 
action between the particles - the interaction is screened. 
At large interparticle separations r, the effective pair po- 
tential between particles cff^^ [r) cx exp(— Ar)/r. Screen- 
ing effects due to the neutralising medium are described 
by an inverse screening length A [l| . Such Yukawa poten- 
tials are ubiquitous in charged colloidal systems arising in 
the classic linearised Poisson-Boltzmann or Derjaguin- 
Landau- Verwey-Overbeek (DLVO) theories for the ef- 
fective potential between spherical charged colloids in 
solution 01 ■ Other systems in which effective Yukawa 
potentials arise include those between dust grains in 
plasma systems |^ and some proteins in solution . The 
Yukawa potential arises not only in diverse physical prob- 
lems, but is also interesting in its own right j^. Note 
that the motivation for considering effective potentials in 
complex systems is driven by the fact that a theoreti- 
cal treatment of the full mixture is often very difficult, 
particularly when there is a big size asymmetry between 
different components of the mixture. It is profitable to 
integrate out the degrees of freedom of one or more of 
the components, thereby incorporating formally their in- 
fluence into an effective Hamiltonian for the remaining 
particles 0. 

In this Letter, we consider a simple model for binary 
mixtures of charged particles suspended in a neutralising 
medium. We examine a dilute suspension of Yukawa par- 
ticles in a solvent with correlation length ^ < A~^ and we 
investigate the (additional) screening effect that the sol- 
vent particles have on the effective interaction between 
the solute Yukawa particles. Specifically, we consider a 
class of binary mixtures in which the pair potentials are 
of the form 4>ij(r) = (/'fj(r) -|- exp(— Ar)/Ar, where in- 
dices i,j = 1, 2 label the two different species and where 
(^fj(r) is a contribution to 4>ij{T) that is short-ranged in 
comparison with the Yukawa term. The amplitudes tij 
depend on the charges on the particles. For example, in 
modelling a binary mixture of charged spherical colloids, 
we would take (/>fj(r) to be a hard-sphere interaction po- 
tential and tij cx ZiZj, where Zi is the (renormalised) 



charge on the colloid of species i However, in what 
follows directly we will not make a particular choice for 
(t>ij{r) or for the sign of e^. We denote species 1 as 
the solvent and species 2 as the solute. We focus on 
the limit of the number density of the solute particles 
P2 ^ and investigate the form of the effective poten- 
tial (f>22'^ (r) between the solute particles in this limit. In 
particular, we obtain the general result that when en, 
£12 7^ 0, so that there is a Yukawa tail cx exp(— Ar)/r 
present in <^ii(r) and (t)i2{r), 4>22'^ (r) decays faster than 
4>22{'r) as r — > CX3 provided the solvent bulk correlation 
length ^ < A~^. On the other hand, when en = ei2 = 
and there is no Yukawa tail in 0ii(r) and 012(7"), the 
effective potential (jfj/ {t) — s- 4>22{t) for r — > cx). This 
Yukawa 'screening' effect is independent of the sign of 
the 'charges', i.e. independent of the signs of e^. The 
results can be re-stated in terms of the solvent-mediated 
potential W22{r), defined via 022 ''^(^) = 022 (j") + W22{r) 
6] . W22 (r) depends on the nature of the solvent and on 
the solvent-solute interaction. If W22{r) faster than 
022(''), as r — + 00, then we describe the solute-solute in- 
teraction as 'unscreened'. However, if VF22('') partially or 
completely cancels the bare potential 022('') at large r, 
we describe the solute-solute interaction as 'screened'. In 
certain mixtures where the bare interactions are purely 
repulsive we find 'super-screening', i.e. (f>22'^ {r) can be 
attractive. 

We may determine the effective interaction between 
two solute particles at infinite dilution from the solute- 
solute radial distribution function .922(7') using the well- 
known result 

mi^{r)^-\n[g22ir)], (1) 

where (3 — (ksT)'^^ is the inverse temperature and .922(7") 
is evaluated in the limit p2 0. For large r, when 
(3(1)22^ [r) is small, it follows that 

/30g^(r)~-/j22(r), r^cx. (2) 

where hij{r) = gij{r) — 1, and in order to determine the 
asymptotic behaviour of 022'^('') ^® must ascertain that 



2 



of /i22 (''') • The total pair correlation functions hij (r) are 
related via the Ornstein-Zernike (OZ) integral equations 
to a set of pair direct correlation functions Cy (r) . The 
OZ equations can be solved formally in Fourier space and 
the solution written (for arbitrary concentration) as 



the solution of Di{q) = 1 — picli{q) = 0, with the small- 
est imaginary part ao and the bulk solvent correlation 
length = 0(7^ 0- Consider now the decay of hij{r) in 
the full mixture. Making the separation of Cij{q) given 
by Eq. © and substituting into Eq. we obtain 



K,{q)^N,,{q)/D{q), 



(3) 



where hij (q) denotes the three-dimensional Fourier trans- 
form of hij{r). The three functions share the same de- 
nominator 

D{q) = [1 - PlCll(g)][l - P2C22{q)] - PlP2Cl2{qf, (4) 



where pi is the density of species i. The numerators 
Nij{q) may be obtained from Refs. ^ or 0- Taking 
the inverse Fourier transform of Eq. evaluating the 
integral by a contour integration around a semicircle in 
the upper half of the complex q plane and assuming that 
the singularities of hij (q) are simple poles one can express 
hij (r) as a sum of contributions from the set of poles at 
{qn} in the upper half of the complex plane 0: 



rh^j{r) = ^ A'lj exp(i(j„r), 



(5) 



where Afj is the amplitude associated with the pole at 
qn- The poles are obtained from the set of solutions 
of D{qn) = 0. The amplitude A^j is related to the 
residue R^^ of qN,j{q)/D{q) by A^^ = R%l2i:. The 
poles are either purely imaginary, 2_= io^Oi or occur as 



a conjugate pair q 
inary pole gives a 
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monotonic contribution to rhij (r) 



of the form exp(— ao'")- A conjugate pair of poles 
gives a damped oscillatory contribution of the form 
2 Aij exp(— aor) cos(air — dij\ where Ai^and Qij denote 
the amplitude and phase respectively In general, 
there are an infinite number of poles. However, the 
asymptotic decay r ^ oo is determined by the pole(s) 
with the smallest imaginary part ao (or do). 

Away from any critical points, the direct pair corre- 
lation functions Cyfr) are known to decay as Cijir) ~ 
—p(j)ij{r), r —f oo When 4>ij{r) has a Yukawa contri- 
bution it is convenient to separate Cij{r) in the following 
way: Cij(r) = cfj{r) — Peij exp{—Xr)/ Xr, which defines 
cfj(r), the short-ranged piece in Cij{r), dependent on the 
form of 4>fj{r) and on the state point. In Fourier space 
it follows that 



c^J(9) =c?'-((z)-a,,/((z2-)-A2), 



(6) 



where a^ = AnPeij/X. 

We seek the pole(s) in hij{q) with smallest imaginary 
part. First, we consider the case when there is no Yukawa 
contribution to the solvent potentials, i.e. an = ai2 = 0. 
The asymptotic decay of /iii(r) in the pure solvent of 
species 1 is determined by the pole(s) in /iii(q), given by 



D{q) ^ A{q) + p2a22Di(q)/p, 



(7) 



where p — q^ + X^ and A{q) = Di{q)[l — /52C22('Z)] — 
PiP2[cY2{q)Y ■ The equation D{q) = has a solution 
q = i[X^ + p2a22Di{q) / A{q)Y/'^ ^ which implies that in 
the limit /92 ~> there is a pure imaginary pole at 
q — qi = iX, provided the ratio Di{q)/ A{q) remains fi- 
nite in this limit. Since we have assumed that there is 
no Yukawa contribution in (j)ii{r), there is no pole at qi 
for the pure solvent and Di{qi) is non-zero and finite. If 
in addition, we assume that both 0^2(91) and C22('3'i) are 
finite it follows that in the limit p2 0, Di{q)/A{q) 1 
for q ~ qi. The amplitude of the contribution to hij (r) 
from this purely imaginary pole at qi is given by 0: 
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(8) 



where the prime denotes the derivative with respect to 
q. It is straightforward to show that in the limit p2 — > 0, 
D'{qi) ~ —2qiDi{qi) / p2a22- Using this result and eval- 
uating the numerators Nij{qi), we find that the ampli- 
tudes of the contributions to hij (r) from the purely imag- 
inary pole at qi = iX are: 
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(9) 



-^+OiP2). 



These results obey the rule AI2 = ^11^22, which general 
considerations demand . Note that in the limit p2 ~> 
the amplitudes ^11,^12 — > 0. Thus the contributions 
from the pole at iX to the decay of /iii(r) and hi2{r) 
are vanishingly small, as one would expect on physi- 
cal grounds. Note further that in the same limit, ^22 
tends to a non-zero constant value — a22/47r = —^£22/^ 
which is independent of any properties of the solvent, 
i.e. the pole at qi = iX gives a contribution to r/i22('') 
of the form — (/3e22/A) exp(— Ar). Thus from Eq. (0) one 
finds that as r — s- 00, 4'22'^ {^) ~^ 022 (t")- In summary, 
when an = ai2 = 0, the effective interaction between 
the solute particles is identical to the bare interaction as 
r 00; the solute-solute interaction is 'unscreened'. Of 
course the argument we have presented supposes that the 
pole qi = iX is the leading-order one; i.e. has the small- 
est imaginary part. If the pure solvent has a correlation 
length ^ > then one expects the asymptotic decay of 
all three correlation functions rhij{r) ~ exp(— r/^) and 
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^'^•i longer ranged than (j)22 {t) . Henceforward we 

restrict consideration to cases where ^ < A~^. 

We turn now to the more reahstic case when an, ai2 ^ 
0, i.e. all three pair potentials have Yukawa tails. Pro- 
ceeding in a similar manner as above (see also the Ap- 
pendix in Ref. the denominator function I^J takes the 
form D{q) = a + b/p + c/p^, where a = [1 — picf^ — 
P2ci5('?)] - PlP2[cf^(g)]^ 6 = [1 - Picl\{q)]p2a22 + [1 - 
P2ci2('7)]piaii + 2/9i/92cf2(g)ai2 and c = pip2(aiia22 - 
0^2)- One set of solutions to the equation D{q) = is 
given by p± = —{b±^/b'^ — 4ac)/2a. This leads to purely 
imaginary poles at q± = ia^ = i-\/A^"^-p±, provided we 
assume that the functions cfj(g) are well-behaved (finite 
and differentiable) on the imaginary axis around q± . The 
leading order pole corresponds to p_ , and in the limit of 
vanishing density p2 — (i.e. c — 0) there is a pole at 
q = iX, so that rhij{r) ^ exp(— Ar), r —^ 00. For 
small concentrations of species 2 we Taylor expand p- in 
powers of c and find the leading order pole is given by 



0.2 



"° " 2X^1 ^2 + 0(p2)- (10) 

Using Eq. ||HJ) we calculate the amplitudes A~j of the 
contributions from this pole to the correlation functions 
/iy(r): 
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FIG. 1: The effective solute-solute potential 4'22^ (r) for sol- 
vent density piX~^ = 3 and P2 — > calculated using the HNC 
closure. Solid line: system A, the point Yukawa mixture, 
with S — and dot-dashed line with non-ideality parameter 
S = 0.1. The dashed line is for system B (5 = 0) where the sol- 
vent potentials decay faster than for the Yukawa case. Results 
are compared to the bare potential <l)22{r) (dotted line). The 
inset displays the radial distribution functions gii{r) (solid 
line) and (712 (r) (dashed line) for 5 = 0. On this scale, there 
is no visible difference between results in systems A and B. 
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Note that the coefhcients of the leading order terms are 
independent of c|J(r) and that the amplitudes obey the 
rule {^12)'^ — ^n^22- 1^ the simplest model of a 
charged system one expects 0^2 — ol\\ol22 since a,,, cx 
ZiZj, the product of the charges on each species. We 
refer to this situation as the 'ideal' mixing rule 8]. Then 
the coefficients of the leading order terms in Eq. Hll|) 
vanish identically and for all three /ly (r) the amplitudes 
corresponding to the pole at q\ = iX will be zero in the 
limit p2 0. Thus the asymptotic decay of /i22(?') (and 
therefore of 022''^ (r)) is determined by the next order pole, 
which generally has a^^ = £, < A~^, so that r4>^2'^ {r) de- 
cays as exp(— r/^), i.e. /oster than 022(''), so that the in- 
teraction between species 2 solute particles is 'screened'. 
In the physical systems where 0^2 7^ 0^110^22, there will 
still be partial screening since one expects the difference 
aiick22) to be small and then (1)^2^ {r) will decay, 
* 00, with the same exponential decay length A~^ 
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as the bare potential, but with a reduced amplitude pro- 
portional to (q!^2 ~ Q^iiQ!22)- We now display results for 
a number of model systems that confirm these general 
predictions. 



The first class of systems is a mixture in which the 
particles interact via purely repulsive (point) Yukawa 
pair potentials for which </'|J(r) = and (pijir) — 
dj exp{—Xr)/Xr 0, 0|. We denote this system A. The 
structure and phase behaviour is described in Ref. Q . In 
this system the effective potential 4>22'^ (r) decays faster 
than the bare potential 022 (^) refiecting the Yukawa de- 
cay of the solvent (species 1) potentials - see Fig.^ We 
compare with system B in which 4'22{f) is the same as 
in system A but where the solvent potentials are mod- 
ified slightly to hasten the Yukawa decay at large r, 
i.e. 0ii(r) — eii exp(— Ar — 7(Ar)-'^")/Ar and 4'i2{r) — 
ei2exp(— Ar — 7(Ar)^*')/Ar. In both systems we choose 
pair potential parameters (3tii ~ 1, /3e22 = 4 and the 
mixing rule ei2 = (1 + (5)^/eiie22 where the parame- 
ter 5 measures the degree of non-ideality 0. The state 
point has solvent density piX~^ = 3 and solute density 
P2X^^ = 10^^, corresponding to the dilute limit. In sys- 
tem B we choose 7 — 10^^, which is sufficiently small 
that the solvent radial distribution functions .911(7") and 
(712(7]) are almost indistinguishable from those in system 
A [ifll- Fig. n displays (/)22'^(^) calculated using the HNC 
closure to the OZ equations 1], which is expected to 
be a very reliable approximation for this model fiuid |^ . 
For system A with 6 = (solid line) there is 'screening'; 
<j>22^ (r) is much shorter ranged than the bare potential 
022(7") (dotted line). In system B with 6 = (dashed 
line) (j>22'^ (r) is indistinguishable from that in system A 
at small r, but for Ar > 2 the results differ significantly 
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FIG. 2: The effective potential (^22 (^) fo'" solvent den- 
sity p\a^ — 0.5 and p2 — > calculated using the PY clo- 
sure for hard-core Yukawa mixtures. Dashed line: system C 
(/3eii = 0.1, /3e22 = 0.4, /3ei2 = 0.2), dot-dashed line: system 
D (/3ei2 = —0.2) and solid line: system E (/3eii = 0, /3ei2 = 0, 
/3e22 = 0.4). In systems C and D, (jjii^i^) is shorter ranged 
than the bare potential 022 (r) (dotted line) but in system E, 
4>22^i''') — <t'22{r) for r > 4a; there is no screening. The inset 
displays the corresponding total correlation function /iii(r) 
for systems C, D and E. 

- there is a maximum near Xr = 3.3 and for Ar > 4, 
4'22^ {^) ~^ 022 (^); the solute-solute interaction is 'un- 
screened'. For S = 0.1 in system A (dot-dashed line) 
4>22^ has a pronounced minimum near Ar — 0.9 and 
is attractive for larger r as the amplitude > - see 
Eq. l(TT|l . We refer to this scenario as 'super-screening'. 
It is remarkable that such effective attraction arises in a 
system where all the bare interactions are purely repul- 
sive and this constitutes a dramatic signal that for S > 
the fluid exhibits liquid-liquid phase separation when the 
densities of the two components are sufficiently high 0. 

The second class of systems is that in which the par- 
ticles have a hard core interaction mimicking (charged) 
spherical colloids. In system C the pair potentials are 
of the form (pijir) — (f>Hs{f) + Eij exp(— Ar)/Ar, where 
4>Hs{r) ~ 00 for r < a and for r > a, the hard-sphere 
diameter jl]|. We choose the parameters A = 0.3cr~"'^, 
/3eii = 0.1, /3ei2 = 0.2, (3e22 = 0.4, corresponding to the 
case where the sign of the 'charge' is the same on both 
species but the magnitude Z2 = 1Z\. We compare with 
system D in which the 'charges' on the particles have 
the same magnitude as in system C but the opposite 
sign, (5t\2 ~ —0.2; all the other parameters remain the 
same. The final comparison is with system E in which 
the solvent (species 1) particles are neutral hard spheres: 
£11 = £12 = while (^22 (t) the same as in systems C and 
D (with /3e22 = 0.4). In Fig. [3 we display the results for 
^'^.^^ ir) calculated using the Percus-Yevick (PY) closure 
to the OZ equations 1] for these three systems at the sol- 



vent density pxa^ = 0.5 and solute density p2(T^ = 10~^. 
PY should be reasonably accurate for such hard core sys- 
tems with additional weak Yukawa tails. We see that in 
system C (dashed line) and D (dot-dashed line) (p^i'^ (r) 
decays much more rapidly than in system E (solid line) 
where (f>22'^ (r) ~ <t'22{i') for r > Aa, i.e. the solute-solute 
interaction is unscreened. 

'Screening' occurs in both systems C and D. Perhaps 
surprisingly, the screening effect is stronger in system C 
where the solvent particles have the same sign 'charge' as 
the solute particles. Note that in case D where the solute 
particles have the opposite sign 'charge' to the solvent 
particles, 4>22'^ {r) can become quite large and positive at 
small and intermediate values of r, as the magnitudes of 
the charges are increased. In this case the screening effect 
discussed here, which applies to the ultimate r ^ 00 
decay of (f>22'^ (r), becomes visible only at very large values 
of r. 

In this Letter we have shown that in Yukawa mixtures 
the interaction between particles of one species can be 
'screened' by the other species. The screening effect does 
not depend on the sign of the 'charges' on the particles 
- see Eq. (|11|) - one observes 'like charge screening', in 
contrast with standard Coulombic screening. Our results 
may have implications for the effective interactions in 
models of charged binary colloidal suspensions. 
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